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' Abstract 

\^ , We calculate the impact factor coupling a virtual photon to a gluon via a massive quark-antiquark 

pair at LL order, but with the imposition of the correct gluon kinematics. Exact analytical results are 
| presented in triple Mellin space with respect to scaled Bjorken x, gluon transverse momentum and heavy 

. quark mass. The application of these results to the calculation of approximate NLL coefficient functions 

needed to relate structure functions to the BFKL gluon is presented. The NLL effects with running 
coupling are seen to lead to a suppression of the small x divergence when compared with the fixed and 
1 running coupling LL results, but less than in the massless case. 

a: 

1 Introduction 

Two factors affecting the successful description of structure function data in deep inelastic scattering are 
the production of heavy quark flavours in the final state, and high energy corrections to the perturbation 
expansion due to the BFKL equation . A recent global fit with LL order small x resummations including 
running coupling effects has shown evidence that resummations can improve the fit to current small x data, 
but also that one should include higher orders in both the fixed order and resummed QCD expansions. For 
a complete NLL analysis one needs the corresponding impact factors coupling an incoming virtual photon 
to the BFKL gluon via a quark pair. Even the massless results are at present unknown, although work is in 
progress The massive results will certainly not be available in the near future. However, it is possible 

to estimate the missing NLL information by supplementing the LL results with information describing the 
correct kinematic behaviour of the gluon. The LL massless impact factors with exact gluon kinematics were 
calculated in [5] and in [S] were shown to approximate well the missing NLL information. The purpose of this 
letter is to extend the calculation of jS] to the massive quark case. This is of phenomenological importance 
as it allows an estimate of NLL effects in the heavy gluon coefficient functions which relate the BFKL gluon 
to measured structure functions. Furthermore, the result has a bearing on the interpretation of the dipolc 
model with heavy quarks beyond LL order (see |7j for a discussion). 
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2 Calculation of the Impact Factors 



The presentation closely follows that of the massless case presented in 5 . Hence we list only the main 
steps here and omit some details, referring the reader to the previous paper. The total virtual photon-gluon 
cross-section is given in 8 . Taking care to separate the transverse and longitudinal contributions correctly, 
one has: 

a L , T ee ^F L , T = J ^T N J ^(7)^t(^V,7,Q 2 /M 2 )(Q 2 ) 7 - 1 , (1) 
where 3(7) is a double Mellin transform of the gluon density, given by: 

«("•■*'» = S 's^w < 2 > 

and <tl,t is a function of the scaled Bjorken x variable: 

x' = x(l + (3) 



with kinematic limits < x 1 < 1, which takes into account the fact that the photon-gluon centre of mass 
energy must exceed the threshold for heavy quark pair production. Hence, we take moments with respect 
to x' rather than x. The variable x g appearing in equation J5J is the longitudinal momentum fraction of the 
gluon that couples to the photon, at the top of the BFKL ladder. This is not the same as the Bjorken x 
variable, which corresponds to the bottom of the ladder. One may relate them using 8 : 

„_ I «P-ft-.)»r+»+fr+"i i #_,(,_,„., (4) 



From these definitions we obtain 



AM 2 \ N f 1 r /i,2\h-2) 

S L (N, 1 ,Q 2 /M 2 )^4aa s [l + -^ r j J dz[z(l z)] 2 J dk 2 ^~ 

d 2 p ( ^ \ K N - (5) 

\p 2 + Q 2 +M 2 {p-k) 2 + Q 2 + M 2 J 




St(N,j,Q 2 /M 2 ) = 

IP 

The kinematic factor 



in equations (1 ~> I ( i fl is the same as in |S] up to the simple addition of M 2 in the denominator. This factor came 
about in JS] due to the relationship between x and x g and the fact that the Mellin transformation of the 
cross-section is taken with respect to the former while that of the gluon is taken with respect to the latter. 
However, for the case of massive quarks we have the additional factor of 



4M 2 ^ N 



Q 2 



(8) 



2 



due to the use of x' rather than x as the real conjugate variable to N. Hence, we have an overall effective 
kinematic factor 

(q 2 + 4M 2 ) N 



\(p - (1 - z)k) 2 + Q 2 + k 2 + M 2 ] N 
representing the difference between the cross-section variable x' and the gluon momentum fraction x g . 



(9) 



One may rewrite equation (J5J as: 



S L (N n ,Q 2 /M 2 ) = 8aa s (Q 2 ) 2 ^ 1 + 



4M" 



where: 



a = J\z[z(i-z)] 2 J ^(k 2 rjd- 



1 



(A-B) 



-K 



D 



dz[z{l - z)f 



(k 2 r / d 2 P - 



(p 2 + Q 2 + M 2 ) 2 ' 
1 



N- 



(p 2 + Q 2 + M 2 )[{p - k) 2 + Q 2 + M 2 ] 
The integrals over k and p can be evaluated using the procedure of [5]. The results are: 



-K 



N ■ 



(10) 

(11) 

(12) 



T(N)(N--/ + 2) 



*(!-*)+ 



7-W+2 



JO 

2 F 1 ( 7 -l,7V-7 + l;^-7 + 3;z); 

JV+2 



NT(N)(N- 1 + 2) W ' J 17 

7-iV-2 

2*1(7- l,JV-7 + 2;iV-7 + 3;z) 



z(l-z) + 7* 



M 2 



(13) 



(14) 



where 2Fi(a,b; c; z) is the hypergeometric function. It does not seem possible to perform analytically the 
remaining integration over the (anti-)quark momentum fraction z as in the massless case. 



Similar formulae are obtainable for the transverse impact factor. First write: 



S T = 2aa s [{DT) + (DT)' - (NDT) - (NDT)'}{Q 2 ) 1 -~ f ( 1 



4M S 



where: 



(DT) = J dzA T (z) J dk 2 (k 2 r- 2 J d 2 p 
(NDT)' = f dzA T (z) [ dk 2 (k 2 y- 2 f d 2 p 



-K 



N : 



(p 2 + Q 2 + M 2 ) 2 

p-(p-k) 



(DT)' = M 2 / dz / dk 2 (k 2 y- 2 / d 2 p 



(p 2 + Q 2 + M 2 )[(p - k) 2 +Q 2 + M 2 1 
1 



-K 



jv. 



(p 2 + Q 2 + M 2 ) 



-K 



N 



(NDT)' = M 2 / dz / dk 2 (k 2 f- 2 / d 2 p- 



p 2 + Q 2 + M 2 )[(p - k) 2 + Q 2 + M 2 ] 



-K 



N 



(15) 

(16) 
(17) 
(18) 
(19) 
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and At{z) = z 2 + (1 — z) 2 . Now introduce also the integrals: 



-K 



N- 



J 

L 



(p 2 + Q 2 + M 2 ) 2 

Q 2 + M 2 



"fc4 



(k 2 r / d 2 P 



[p 2 + Q 2 + M 2 } [(p - k) 2 + Q 2 + M 2 } 



(p 2 + Q 2 + M 2 )[(p - k) 2 + Q 2 + M 2 ] ' 



K N ; 



-K 



N ■ 



so that: 



(DT) - (NDT) = [ dzA T (z) 
Jo 



J(z)-I(z) + -L(z) 



(20) 

(21) 

(22) 
(23) 

(24) 



Each of the integrals in equations l|18l22fl can be calculated by analogy with the longitudinal case. The 
results are: 



mrv - (o 2 \^ M2 ^r( 7 -i)r(JV-7 + i) 

{ ' ~ W ' Q 2 T(N)(N-y + 2) 



z(l -z) + 



M 2 



2-1-T-JV-2 



z N -^+\l - z) N 2^(7 - 1, N - 7 + 1; N - 7 + 3; z) 



{NDT)' = (Q 2 y 



_ x M 2 27rr(7 — i)r(iv -7 + 2) f 1 



dz 



z(l-z) + 



M 2M-N-2 



Q 2 r(7V + i)(iV-7 + 2) J 
x ^-t+ 2 (1 - z) N 2^1(7 - 1, W - 7 + 2; JV - 7 + 3; z) 



(25) 



(26) 



dzA T (z)I{z) = {Q 2 Y 



. 1 7rr(7-l)r(JV-7 + l) 
T(N){N- 1 + 2) J Q 



dz 



z(l-z) + 



M 2 

Q 2 



2-17-jv-i 



z JV- 7 +i (1 _ z) jv tl _ 2z(1 _ z )] 2 Fi(7 — 1,JV — 7+l;iV — 7 + 3; z); (27) 



/' 

Jo 



dzA T (z)J(z) = (Q 2 )^ 



_ 1 2 7 rr(7-l)r(jV-7 + 2) 
T(iV+ l)(JV-7 + 2) 



rfz 



,1 N M 



7-AT-l 



z jv- 7 +2 (1 _ _ 2z (i — z)] 2^1(7 — 1,AT — 7 + 2;iV — 7 + 3; z); (28) 



dzA T (z)i(z) = (Q 2 )T- 



27rr( 7 )r(^- 7 + l) f 1 



r(7V + l)(iV-7 + l) 7o 



dz 



z(l-z) + 



M 2 



7-iV-l 



z N -i +1 (l - z) N [l - 2z(l - z)] 2 Fi( 7 , JV - 7 + 1; JV - 7 + 2; z). 



(29) 
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Putting things together, the total impact factors may be written as: 



<? (N -v C*IM»\ 87raa s r(7-l)r(JV-7 + l) / 
S L (i\T, 7 ,Q /M ) = T{N){N _ 1 + 2) I 1 



4/VP 



2\ N /•! 



*(1 - *) + w 



7-JV-2 



2 (TV — -Y + 11 2 

2-Fi(7-l,iV-7 + l;iV-7 + 3;*)--i J 2-^1(7 — 1, TV — 7 + 2; TV — 7 + 3;z) 



7, Q /M ) r(A0(A r_ 7 + 2) ^ + -Q2 

M 
Q 

2(TV - 7 + l)z 



N 

2\ N /"I 



*(1 - Z) + 



7-JV-l 



(30) 
( 1 - z) 1 



X { z(l- z) 



2-[z(l-z) + „ 2 



N 



2 F 1 (<y-l,N-~/ + 2;N -7 + 3;*) 



2*1(7 - 1, TV - 7 + 1; N - 7 + 3; z) 
( 7 -l)(JV-7 + 2) 



iV(iV-7 + l) 



[l-2z(l-,z)] 2J Fi(7,iV-7 + l;JV-7 + 2;«) 



(31) 



We also note that one can use equation (jSJ) to rederive the known LL impact factor in a simpler form. 
Setting TV = in this equation, one may carry out the integrals over k 2 and p as in the massless case 
presented in ^J] to give: 



S L (0, 7, Q 2 /M 2 ) = Anaa s ^—^ljL^ f dz[z[l _ z) f 



r(2 - 2 7 ) (3 - 2 7 ) J 
The integral over z can be performed after making the change of variable: 



7-2 



(32) 



(1 - yT=u )/2, 2<l/2 
(l + VT _ ^)/2, 2>l/2, 



followed by u' = 1 — u to get: 



£ L (0,7,07*0 



27raa s T 3 (l - 7 )r( 7 ) (1 - 7) ^ M 
15 r(2 - 2 7 ) (3 - 2 7 ) 



2\ 7-2 



2 \ 7-2 



.1 



4Af 2 / 



*i 2- 7 ,-;-; 



1 7 Q 2 



2 ' 2 ' Q 2 + 4M 2 
(33) 

This can be shown to be equivalent to the previously published result for the massive LL longitudinal 
impact factor 4 [S], which has a slightly more complicated analytic form instead of being proportional to a 
single hypergeometric term. No particular simplification is possible for the transverse impact factor, however, 
as this has more terms in the integral definition ©. 



3 Results in Triple Mellin Space 

It seems that the integrals in the preceding section cannot be performed analytically due to the factors of 
[z(l — z) + M 2 /Q 2 ] in the integrand. However, one may decouple the mass dependence and obtain closed 
analytic forms for the impact factors by taking a further Mellin transform with respect to M 2 /Q 2 , defined 
by: 

^,MW)=(l + iV. 7l ). (34) 

4 Note that a factor of 7/(47r 2 o) is needed to convert to the notation of the impact factor /i_l( 7 , M 2 /Q 2 ) as defined by Catani, 
Ciafaloni and Hautmann. 
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where the kinematic factor of (1 + 4M 2 /Q 2 ) from x' has also been included. Using the result: 



dyy^Ml -z) + y] b = ^i-W-Ti b) _ (3g) 



n-b) 

one finds from equation (|30(l : 

5 , 87raa s r( 7 - l)r( 7l - l)r(JV - 7 + l)r(JV + 3 - 7l - 7 ) f 1 



dz 2^(1 -^)T+7i-l 



r(7V)r(7v - 7 + 3) J0 

s 2(AT-7+l)z 

2*1(7 - 1, JV - 7 + 1; JV - 7 + 3; z) - -i ^ 2*1(7 - 1, N - 7 + 2; JV - 7 + 3; z) 



iV 



(36) 



The integral over z now reduces to a conventional hypergeometric integral, and one may use the standard 
result [Trl] : 

dxx a - l (l - xf- 1 2 F 1 (a,b;c;x) = 3*2(0, 6, a; c, a + f3; 1) 



to yield: 



r(a + /9) 

87raa5r(iV + 3 — 71 — 7)r( 7l — l)r(7 + 71) 
f(A0 

3 G 2 (7 - 1, N - 7 + 1, 71 + 1; -/V - 7 + 3, 7 + 2 7i + 1; 1) 



(37) 



N 



3 G 2 ( 7 ~ 1, N ~ 7 + 2, 71 + 2; TV - 7 + 3, 7 + 2 7l + 2; 1) 



(38) 



where we have used the Meijer G function 36*2(0, 6, c; d, e; z) = r(a)r(b)r(c)/[r(d)r(e)] 3*2(0, 6, c; d, e; 2). 
One can check this result against the massless result for the impact factor after the integration over z has 
been performed 5 . The limit M 2 — > corresponds to 71 — > 1 according to the Mellin variable definition of 
equation and so one has: 

dM 2 /AA 2 \ 7l ~ 2 
lm ' ' "qHq 2 ") ^(^7,M 2 /Q 2 )=^(7V,7,7i)l 7l ^i 



dM 2 /Af 2x71-2 



lim 



[r( 7l -i)s L (iv,7,o)] 



71— >1 • 



Thus equation (|3T))l . divided by T(7i — 1), reduces as 71 — * 1 to the corresponding result in the massless 
case [5]. When M — a simplification occurs, and one can eliminate the Meijer G functions in favour of ip 
functions. Such a simplification is not possible in the massive case of equation (|38|1 because of the non-zero 71. 
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Similarly, one can Mellin transform equation 1|31[) to 71-space, with the result: 
2naa s T(2 - 71 - 7 + A0 r (7i - l)r( 7l + >y)(N - 7 + 1) 



T(N) 

x |2 3 G 2 (7 - 1, N - 7 + 1, 71 + 1; N - 7 + 3, 7 + 2 7i + 1; 1) 
- -^r 3 G 2 (7 - 1, ^ - 7 + 2, 71 + 2; iV - 7 + 3, 7 + 2 7l + 2; 1) 



1 



iV(JV-7+l) 



' 3 G 2 (7, # - 7 + l,7i; JV - 7 + 2,7 + 271 - 1; 1) 



71 + 7 - 1 



2 3 G 2 (7, N - 7 + 1,71 + 1; N - 7 + 2,7 + 2 7l + 1; 1) 



(2 - 71 - 7 + AT) 



2 

— 3 G 2 (7 - 1, iV - 7 + 2, 71 + 1; 7 + 2 7 i, iV - 7 + 3; 1) 



( 7 i+ 7 -l)(iV-7+l) 
- 3 G 2 ( 7 -l,iV-7 + l, 7157 + 271 -l,-/V-7 + 3;l) }■ (39) 

4 Implications for Structure Functions 

The impact factors discussed in this paper contain some of the information in the unknown NLL heavy impact 
factors. In the massless case, we have found evidence that the exact kinematics impact factors approximate 
well the true NLL calculation [S]. The comparison relied upon knowledge of the fixed order coefficient and 
splitting functions P qg and C^g up to NNLO ^2 ^] . A similar analysis is not possible for the heavy impact 
factors, as the NNLO heavy coefficient and splitting functions are not known, and thus there is no subleading 
small x information at fixed order available for comparison with resummed results. Nevertheless, by analogy 
with the massless case it is a reasonable assumption that the NLL impact factors are approximated well 
by the exact kinematics results and we can therefore use them to investigate the relevant phenomenological 
quantities. First one interprets the impact factors in terms of the coefficient functions relating the BFKL 
gluon density to the proton structure functions using the kx factorisation formula |131 in double Mellin 
space: 

Ff = (7^(7, N, M 2 /Q 2 )g( 7 , AO, (40) 

where Ff 1 (i E {2,L}) is the heavy flavour contribution to the structure function 5 . Here Q(j,N) is the 
unintegrated gluon density, related to the integrated gluon in Mellin space by Q(-f,N) = 75(7, N). Taking 
into account also the normalisation factor in equation (Q, one identifies: 

C L , S ( 7 , N, M 2 /Q 2 ) = -r^-Std, N, M 2 /Q 2 ); (41) 
C 2 , g ( 1: N,M 2 /Q 2 ) = -t^{S T ( 1: N,M 2 /Q 2 ) + S L ( 1: N : M 2 /Q 2 )}. (42) 

Expressions for these quantities in the physical space of x' and Q 2 are obtained by solving the BFKL equa- 
tion. We here adopt the solution method of ^B] to investigate the behaviour of the NLL resummed coefficient 
functions. The calculation demonstrates that the results of equations l|3UI I31fl are of phenomenological use 
even though the final integration in z cannot be performed analytically. 

For a NLL solution of the BFKL equation, one requires the impact factors truncated at O(N) and expanded 
as a power series in 7. Expansion of the integrands of equations (|30l I31|) is possible using the following 

5 Note that this contribution corresponds to a fixed flavour description of the structure function, where the number of active 
quark flavours is the same at all Q 2 . In practice one should use a variable flavour scheme at high Q 2 but we do not consider 
such a scheme here. 
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integral representation of the hypergeometric function 



F 1 (a,b;c;z) = — ffi , / dvv b -\l - v)^ 6 ^! - v Z y a . (43) 



2 -T(bnc-b)J 
Then one has, for example: 

2 F 1 (y-l,N- 1 + l;N- 1 + 3;z) = ^£Lzl±V f dv (1 - «)(1 - vz)\\ + Nlogv + 0(N 2 ) 

T(N -7+1) Jo { 

- 7[log(l - vz) + logv + (logwlog(l - vz) + log 2 v)N + 0(N 2 )] + C(7 2 )|. (44) 
This is a power series in N and 7, of the form: 

2 F 1 ( 1 -l,N- 1 + l-N- 1 + 3;z) = Y,Y. Nn ^ m [ dvf nm (v,z), (45) 

where the integrands fnm(t, z) become progressively more complicated as n and m increase. A few low order 
terms can be integrated analytically. When this is not possible, one can parameterise the integral using a 
suitable function of z: 



[ dv f nm ( V , Z )~J2 k ^ zr + E fc «™ (! - z ) ^o- - *) 

Jo 



(46) 



for some range of r and s. Values of the integral can be taken for various values of z, and fitted to this 
functional form using a least squares routine. The logarithmic terms are sometimes needed to obtain a good 
fit, but are weighted by (1 — z) noting that an expansion about 7 = of equation 144(1 is finite as z — > 1. 
This is also the case for the other hypergeometric functions encountered in the impact factors. 

One can now expand in N and 7 the z integrands in equations (|3(JI 131(1 . after substituting in the parameterised 
hypergeometric functions. Then one has: 

S,( 7 ,iV,M 2 /Q 2 ) = I dzY, \sl°l(z,M 2 /Q 2 ) + NS^(z,M 2 /Q 2 )} 7" + 0(N 2 ) (47) 
J ° n 1 J 

where the power series in N is truncated at NLL order. For given values of M 2 /Q 2 , the z integrals can be 
calculated numerically to obtain the expanded massive impact factor at a given momentum scale. For exam- 
ple, choosing M = 1.5GeV and A = 150MeV for the QCD scale parameter, one finds at t = log(Q 2 /A 2 ) = 7, 
i.e. Q 2 = 25GeV 2 : 

C Lg (7, N, M 2 /Q 2 )\ t=7 = ^[.5184 + .2069 7 + .9468 7 2 + .3799 7 3 + 1.155 7 4 + A(-.4869 - 1.105 7 - 1.936 7 2 

47T 

- 2.880 7 3 - 3.7857 4 )]; (48) 
C 2g (7, N, M 2 /Q 2 )\ t=7 = ^[3.318 + 1.280 7 + 7.095 7 2 + 3.473 7 3 + 10.35 7 4 + A(-5.154 - 8.613 7 - 20.22 7 2 

47T 

- 28.41 7 3 - 45.957 4 )]. (49) 

One may check the LL coefficients in equations (1481 149(1 agree with the results of expanding the LL impact 
factor found in Furthermore, one may verify the NLL part of Cl b for M 2 /Q 2 — > against the 

massless exact kinematics result 0. One cannot do this for C 2 ,g due to the fact that the transverse impact 
factor diverges as M 2 /Q 2 — > 0. However, it is possible to ascertain the correct asymptotic behaviour. At LL 
order one has: 

Q^(7,M 2 /Q 2 )|^ — ► A( 7 ) J +/ 2 ( 7 ). (50) 
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At NLL level and beyond, one can make a similar ansatz by promoting /, = ^(7, TV). Lesser powers of 
M 2 /Q 2 would vanish as M 2 /Q 2 — ► 0, and higher powers are inconsistent with the fact that the fixed flavour 
coefficient functions in Q 2 space contain at most collinear divergences ~ a 7 g log"(M 2 /Q 2 ). From equation 
1(50)1 one finds: 



1 



C 2 , s ( 7 ,iV,M 2 /Q 2 ) + 



M 1 



dC 2 ^(j,N,M 2 /Q 2 ) 



d log Q 2 /M 2 



A (7, N) = Jim [C 2 , g (rr, N, M 2 /Q 2 ) - / 2 ( 7 , N)] 



There is also the consistency relation: 

1 C 2 ^ g ( 1 ,N,M 2 /Q 2 ) 



dC 2 , g {i,N,M 2 /Q 2 
d\og{Q 2 /M 2 ) 



lS' 2 ( 7 ,N), 



(51) 
(52) 

(53) 



where £"2(7, iV) is the massless impact factor with exact gluon kinematics. Thus one has / 2 (7, A^) — 
S'2(l,N) as an all orders result in 7. Then f 1 can be found as a series expansion from equation (|52fl by 
substituting the power series for C 2 . g (7, N, M 2 /Q 2 ) at some sufficiently low value of M 2 /Q 2 . It can be 
verified against the result obtained in the LL case. Not only is this an important check of the calculation, 
but the asymptotic limit is of phenomenological importance for disentangling the heavy flavour coefficients 
in a variable flavour scheme at high Q 2 (see for how this works at LL order). 

Armed with the series expansions of equations ((481 149(1 , one can use the method of JHj to obtain the cor- 
responding estimated NLL coefficients Ci t9 in x' and Q 2 space. Results are shown in figure 01 6 . We also 
show the LL results with running coupling corrections and the strictly LL results. Note that all of these 
coefficients are in the Qo scheme, where the gluon is defined by the solution of the BFKL equation. In 
principle a transformation is needed to obtain the corresponding results in the more commonly used MS 
scheme. Looking at the figure, one can see that the LL coefficients are strongly divergent at small x. This 
divergence is tempered significantly by the inclusion of the running coupling (as already noted in ^3|2]) and 
the effect of the approximate NLL corrections is to suppress the small x divergence yet further. A similar 
suppression is observed in the massless case ^SJ, and in [2] it was noted that a softening of resummation 
effects in the moderate x region is needed to achieve a good fit to scattering data. 

We note that a smaller suppression of the small x divergence occurs in the massive impact factor results 
than in the corresponding massless quantities when exact kinematics are included. One can examine this by 
comparing the factor Kfj of equation © with the equivalent factor in the massless case, which is: 



K 



M=0 
N 



(Q ) 



2\N 



[( P - (1 - Z)k) 2 + Q 2 + k 2 } 



2]N 



(54) 



This has a maximum value of unity at k 2 — and (p — (1 — z)k) 2 = 0. It then acts to suppress the integrand 
for the impact factor as k and p move away from these values, becoming noticeable at a scale given by 
(p — (1 — z)k) 2 + k 2 ~ Q 2 . Similarly, has a maximum value of unity at the same values of k and p, 
and z = 1/2. But the typical scale at which a significant suppression occurs is increased by the presence 
of the heavy quark mass to (p— (1 — z)k) 2 + k 2 ~ Q 2 + 4M 2 . Thus there is less of an effect in the massive case. 

This reduction of the suppression due to exact kinematics in the heavy flavour structure functions is of 
particular phenomenological interest. NLO QCD fits have a tendency to undershoot data on F| at low x 
(e.g. mi])- The LL resummed fit of |2] also underestimated F| at very low x, and thus one would hope that 
in a NLL resummation the heavy coefficients are suppressed less with respect to the massless results. It is 
encouraging that this indeed seems to be the case. 

6 There are also corrections ~ 0(A 2 /Q 2 ) to these results due to ambiguities in the derivation 1151 . We assume, however, that 
they are small at t = 7. 
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Figure 1: Estimates for the NLL resummed fixed flavour heavy gluon coefficient functions for Fl and F2 
(solid) using the exact kinematic impact factors at t = 7, rif = 3. Also shown is the LL result with running 
coupling corrections (dashed) and the strictly LL result (dotted). 
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5 Conclusions 



In this letter we have calculated the virtual photon-gluon impact factor with massive quarks and exact gluon 
kinematics, presented in equations (|3UI I31|) . Closed analytical results do not seem to be possible in M 2 /Q 2 
space (there is a remaining integration over the heavy quark momentum fraction z), although are possible 
in the scaled triple Mellin space of section and are given in equations and Furthermore, one 

can still use the results in M 2 /Q 2 space for phenomenology. As an example, we presented estimates for the 
NLL resummed fixed flavour heavy flavour gluon coefficients Ci^ g and The estimated NLL effects are 

seen to suppress the low x divergence when compared with both the fixed and running coupling LL results. 
Such an effect is encouraging, as this is hopefully what is needed to achieve a good fit to structure function 
data over the complete x range. Furthermore, the massive results are suppressed less than the corresponding 
massless results as a consequence of the kinematic constraint on the partonic centre of mass energy at the 
heavy quark vertex. This is consistent with previous global fits at NLO and LL orders in the QCD expansion 
that underestimate the charm data at very low x. Work on implementing our results in an approximate NLL 
global fit is ongoing. 
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